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Abstract 

Deep inelastic scattering is considered in a statistical model of the nucleon. This 
incorporates certain features which are absent in the standard parton model such 
as quantum statistical correlations which play a role in the propagation of particles 
when considering Feynman diagrams containing internal lines. 

The inclusion of the 0{as) corrections in our numerical calculations allows a 
good fit to the data for x > 0.25. The fit corresponds to values of temperature 
and chemical potential of approximately T = 0.067 GeV and fi = 0.133 GeV. 
The latter values of parameters, however, give rise, for all x, to a large value for 
R = aL/cTT- 



UCT-TP 187/92 
BI-TP 92-45 
February 2008 



1. Introduction 



In this paper we further develop a statistical model for the structure functions 
of the nuclcon [1,2,3,4] based on the MIT-bag model [5]. The scattering of high 
energy leptons from a thermalized gas of quarks and gluons confined to the volume 
of the nucleon was calculated by considering lowest-order diagrams in perturbative 
QCD. In refs. [6,7,8,9] these considerations were extended to the first order in the 
strong coupling constant ag since the effect of gluons is absent in zeroth order. 

The phenomenological success of the parton model indicates that the number 
of partons in the nucleon is very large due to the 1/x behaviour of the structure 
functions at small x. 

The parton model neglects quantum statistical correlations due to the pres- 
ence of identical quarks and gluons in the initial and final states. In the statistical 
model this is taken into account through the use of Fermi-Dirac and Bose-Einstein 
distributions for quarks and gluons, respectively. Stimulated emission factors for 
final-state gluons and Pauli-blocking factors for final-state quarks are incorporated. 
The propagation of particles through a many-body medium is taken into account 
by using thermal Feynman rules for propagators and vertices. These effects are 
negligible in lowest order, however, they are important for higher order diagrams 
especially in the small x region. 

In ref. [9], all processes contributing to order as (see Fig. 1) to deep inelastic 
scattering of leptons off a heat bath of quarks and gluons, were taken into account 
exactly. All infrared, collinear and ultraviolet divergences cancel in the framework 
provided by the real-time formalism of finite-temperature quantum field theory. 
The final expressions from these analytical calculations are given in the appendix. 

The presence of the medium leads to a dramatic increase in the length of 
the calculation. For this reason many authors [10,11,12,13, 14,15] use special 
kinematics [e.g. in lepton pair production from a thermalized quark gluon plasma 
they would use a virtual photon with zero spatial momentum: q = (goi^ = 0)]. 
In our case we perform calculations for the case of general kinematics [i.e. q = 

(go, 5 7^0)]. 

Results from analytical calculations of dilepton production from a quark-gluon 
plasma to order ag {q^ > for the photon giving rise to the lepton pair), which is 
closely related to the process of deep inelastic scattering of leptons off a heat bath 
of quarks and gluons (g^ < for the photon stemming from the lepton beam), are 
presented in ref. [9] . In Section 4 we present results from numerical calculations of 
these expressions. There it is shown to what extent recent deep inelastic scattering 
data can be reproduced by the statistical model. 
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2. Deep inelastic electron-proton scattering 



In our treatment of deep inelastic scattering we consider the scattering of a 
virtual photon with 4-momentum q — {u, 0, 0, q^) off the proton (considered as a 
heat bath of quarks and gluons). To order ag the processes depicted in Fig. 1 
contribute. According to Fig. 1 there can be more than two particles in the initial 
state of a reaction so that the rate of reactions (as opposed to a cross section) is 
the natural quantity to consider when relating the measured structure functions 
on the hadronic level with the results from our calculations of thermal averages of 
the fundamental (virtual) photon-quark interaction processes. 

The deep inelastic cross section for charged leptons can be written in such 
a form that it depends linearly on F2 and non-linearly and very weakly on R 
[16]. Experimentally, F2 and R can be disentangled by measuring deep inelastic 
cross sections at the same {x,Q'^) point but at different values of the incident 
beam energy. In practice, as these cross sections depend only weakly on R, all 
measurements of R have large errors. 

To each of the processes (a) to (f ) in Fig. 1 there correspond the two Feynman 
diagrams in Fig. 2. These two diagrams, as drawn in Fig. 2 with exactly the same 
four-momentum assignments and directions of arrows, are used for each of the 
processes (a) to (f) in Fig. 1. When the energy components of the latter four- 
momenta are all positive the two Feynman diagrams refer to process (a) in Fig. 1 
and the energy of a particle in diagrams (b) to (f) in Fig. 1 is considered to be 
negative when it changes from the left hand to the right hand side or vice versa 
relative to diagram (a) in Fig. 1. 

To each of the diagrams (g) and (h) in Fig. 1 there corresponds the six 
Feynman diagrams in Fig. 3. These six diagrams, as drawn in Fig. 3 with exactly 
the same four-momentum assignments and directions of arrows, are used for each 
of the processes (g) and (h) in Fig. 1. For ko > and kg > we are considering 
the three-particle process (g) involving quarks and for ko < and /cq < we are 
considering the three-particle process (h) involving antiquarks. 

The relation between the measured structure functions on the hadronic level 
and the results from our calculations of thermal averages of the fundamental (vir- 
tual) photon-quark interaction processes is discussed in refs. [2] and [7]. One 
obtains on the hadronic level, as an example, for the 0{as) process of gluon emis- 
sion from a quark in Fig. 1(a): 
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according to the momentum assignments in Fig. 2 and where V is the volume 
of the nuclcon. In eq. (1) A refers to the heUcity of the virtual photon. The 
matrix element for gluon emission from a quark for a given polarization, A, of the 
space- like (g^ < 0) photon (7*) is given by 

Ma = Ml + M2 (2) 

with 

Mi=iZ(fc',s')(-^^7A*)[e^(pF^^''(^ + 9)(^e7.)e^(g)^fc,s) (3) 

M2=iZ(A;',s')(^e7^)e^^^''(fc-p)(-^^77.)[eg(p)]*«(fc,s) (4) 

where we suppress reference to colour factors and fractional charges and where iS^^ 
denotes the 11-component of the fermion propagator in the real-time formalism of 
finite-temperature quantum field theory. 

The equations (3) and (4) provide explicit expressions for the matrix element 
in eq. (1) (with coupling constants, colour factors and fractional charges factored 
out into explicit normalization constants). We take Nc = 3 and include up and 
down quarks in our model. 

Fermion and boson distribution functions can be defined in such a manner 
that they are valid for both positive and negative energies. They are for a fermion 

^^(^^) = ^^^^ 

where 

Xk = \ko\ - Ate(^o) (6) 
where e(A;o) gives the sign of ko and for a boson 

^B(|feo|)= ^^11,^1 (T) 

3. Calculation of the structure functions 

The zeroth order in 0:^ contributions can be calculated in a straightforward 
manner [2]. The calculation of 0{as) contributions require a much larger effort. In 
the calculations of the 0{as) contributions, the phase space and loop integrations 
can be analytically reduced [9] to an integration over the two energy variables K 
and po as shown in the final expressions for four-particle processes given by eqs. 
(Bl) and (B2) and for three-particle processes given by eqs. (CI) and (C2). These 
double integrals cannot be analytically calculated. Their numerical calculation is 
discussed in Section 4.2. 
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We distinguish between the zeroth and first order in cXg contributions by 
writing: 

F2 = Fi'^+F^'\ (8) 

Since we are also interested in the structure function R which is expressible in 
terms of F2 and Fi , we also calculate 



F^=F}'^+F}'l 



Zeroth order in as expressions for F2 and Fi in the statistical model are gi 
ref. [2]: 
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where the sum over q is over quark flavours only and where 

_ Mx jlq 
Zq = ^-Y 

and 
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The 0{a.a) parts of expressions for F^ and F\ in the statistical 



model are: 
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with the expressions in curly brackets given in eqs. (CI), (C2), (Bl) and (B2). 
The have to be evaluated for each flavour of quark, since the chemical potential 
appearing in the statistical factors ^4 and S^w differ for each flavour as we sum 
over the flavours of quarks. As explained in Appendices B and C, the K,po 
integration regions for four- and three-particle processes are as given in Figs. 4 
and 5, respectively, and the integrand for four-particle processes is different in each 
subregion shown in Fig. 4. 

4. Numerical Evaluation 

The numerical calculation of the structure functions F2 and R was performed 
in two parts. The zeroth order is discussed in the next section while the contri- 
butions of 0{as) are discussed in Section 4.2. Results from the total calculation, 
i.e., zeroth order plus 0{as) results are presented in Section 4.3. 



4.1. Zeroth Order 

We chose to flt to a plot of F2 versus x, which is determined by parton 
distributions which reproduce a wide range of experimental data. Such parton 
distributions were published in 1990 by Kwiecinski et al. in ref. [17]. They 
parametrize the distributions at Qq = —q'^ = 4 GeV^ and evolve up in in 
order to test their parton distributions against experimental data. We flt to their 
so-called B_ flt at = 4 QeV^. 

The values for temperature and chemical potential which give the best flt for 
F2 versus x are used to calculate the corresponding plot of R versus x which we 
will compare with experimental data given by Whitlow [18]. In the deep-inelastic 
scattering limit, x,T, jj,q flnite with v going to inflnity, the expression for ^2°^ in 
eq. (10) is dominated by the /o term which can be evaluated analytically. In Ref. 
[2] the latter expression was used to flt the experimental data by Aubert et al. 
[19]. 

In this work we keep the full expression for ^2*^'' as given in eq. (10) and flt 
to F2 as determined by the parton distributions by Kwiecinski et al. in ref. [17] 
(their B_ flt). The values for temperature and chemical potential which produce 

the best flt to the plot of F2 versus x are used to calculate F^^^ in eq. (11) too. 

Then the structure function R which is expressible in terms of F^'^ and F^^^ can 
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be calculated. The chemical potential for the down quarks was kept fixed in terms 
of the chemical potential for the up quarks according to 

(19) 

(as was done in ref. [2]) in order to reproduce the ratio of up to down quarks 
as being 2 to 1. This is approximately the case for a quark gas at very low 
temperatures. 

As opposed to the /o function, the functions /i and /2, which also appear in 
the expression for ) in eq. (10), cannot be evaluated analytically. They were 
evaluated numerically. 

For X > 0.4, the zeroth order in as calculation produced a good fit to the 
= 4GeV^ data of Kwiecinski et al. as shown in Fig. 6. 

According to the following rough estimates, the failure of the zeroth order in 
cts theory of the statistical model to reproduce the low x behaviour of F2 could 
be due to the finite volume of the nucleon. From the mass shell condition on the 
four-momentum k' = k + q of the outgoing quark in the zeroth order in ag (lowest 
order Feynman diagram in Fig. 3), one obtains a lower limit on of the incoming 
quark as follows. From the mass shell condition for massless quarks [{k + q)'^ = 0] 
and the relation x = Q"^/ {2Mv) one obtains 

-1 < cOSer.= ^^^^^= < 1. (20) 

2\k\Vp^ + 2Mi^x ~ 
The first inequality in the latter line produces 

Q2 2Miyx 



\k\ > 



^%(!?(^). (21) 



2 

Another lower limit on is obtainable from the Heisenberg uncertainty principle 
of the general form AxAk > 1. Taking Ax equal to the diameter of the nucleon 
(2 fermi = 10.14 GeV~^), one obtains Ak > (10.14)"^. From this and the value 
of the lower limit on \k\ in eq. (21) one could deduce that our model is only valid 
for X values satisfying 

Mx , X 1 , X 

^ (10-14)-', (22) 

i.e., X > 0.21. 

The zeroth order in as results for the structure function R are shown in Fig. 
7. For X > 0.25 these results are consistent with the = 4GeV^ experimental 
data given in Ref. [18]. 

The zeroth order in as calculation is only able to reproduce the data for values 
of X > 0.4. In Section 4.2 we investigate to what extent the 0{as) corrections are 
able to extend the fit to values of x smaller than 0.4. 
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4.2. The zeroth- plus first order in ag calculation 

In this section we take the 0{as) corrections into account by numericaUy 
calculating the quantities given in eqs. (15) and (16) and adding them to zeroth 
order in results according to eqs. (8) and (9). The zeroth order in quantities 
given in eqs. (10) and (11) are calculated in the manner discussed in Section 4.3. 
Due to the modifications introduced by the 0{as) contributions, we search for 
new values of temperature and chemical potential which will produce a fit to the 
plot of F2 versus x of Kwiecinski et al. given in Fig. 6(a). 

We keep the chemical potential of the down quarks fixed in terms of the chem- 
ical potential of the up quarks according to eq. (19). The values for temperature 
and chemical potential which produce the best fit to the data are used to calculate 
Fi in eq. (9). Then the structure function R which is expressible in terms of F2 
and Fi can be calculated and compared to the experimental data. 

According to the 0{as) expressions, ^2^'* and F^^'^ in eqs. (15) and (16), 
one needs to numerically calculate the double integrals in eqs. (Bl), (B2),(C1) 
and (C2). The K^po integration regions for four- and three-particle processes are 
shown in Figs. 4 and 5, respectively. These integrals are devoid of ultraviolet, 
infrared and collinear singularities as discussed in ref. [9]. However, individual 
terms in the integrands can still become infinite along certain lines in the K,po 
integration regions. This happens when factors in denominators or arguments of 
In functions in some terms become equal to zero along certain lines. 

Examples are singularities which arise along the line along which the quantity 
b is equal to zero. The latter singularity arises in the terms proportional to In 16*41 
and IG4I in eq. (Bl). One of the subregions which contain the line along which b 
is equal to zero is subregion D in Fig. 4. 

Even though singularities from such terms cancel each other, the terms indi- 
vidually give rise to singularities. Therefore one should take care not to approach 
the line along which 6 = too closely in the numerical integration. For the same 
reason one should steer clear of other lines along which other individual terms 
in expressions for three- and four-particel processes are singular. This is accom- 
plished by implementing a trapesium integration method whereby the integrand 
is evaluated at the discrete points of a two-dimensional lattice constructed in such 
a way that the two nearest successive lattice points to any of the lines along which 
singularities arise, are always located an equal distance on both sides of the rel- 
evant line. The value of the integral is the stable value attained in successive 
numerical integrations in which, independently and in steps, the density of lattice 
points is increased and the range of integration is extended (The presence of ther- 
mal distribution factors in terms suppresses contributions at large values of the 
integration variables.). 

4.3. Results 

In Section 4.1 it was mentioned that the zeroth order in as calculation is able 
to reproduce the data of Kwiecinski et al. on F2 versus x for values of a; > 0.4 only. 
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From Fig. 8 it can be seen that, by including 0{as) corrections, the reproduction 
of the data of Kwiecinski et al. could be extended to values of x smaller than 0.4. 
A reasonable fit could be obtained for values of x > 0.25. This is compatible with 
our estimate for the limited range of validity of our theory {x > 0.21 as derived in 
Section 4.1) due to the finite volume of the nucleon. 

For the fixed values of the parameters and the fixed expression for the chemical 
potential of the down quarks in terms of the chemical potential of the up quarks in 
eq. (19), a fit to the data of Kwiecinki et al. is possible for values of temperature 
and chemical potential in the immediate vicinity of T = 0.067 GeV and /Xup = 
0.133 GeV only. 

Figs. 9 and 10 show how the results for F2 versus x from the calculation 
including zeroth order plus 0{as) contributions vary as a function of tempera- 
ture at fixed chemical potential and as a function of chemical potential at fixed 
temperature, respectively. 

From Fig. 8 it can be seen that the 0{as) corrections to F2 is negative 
for large x. Such negative values could in principle arise from interference terms 
from three-particle processes since, by neglecting the 0{a1) contributions to the 
amplitude squared of the 0{a^) and 0{al) Fcynman diagrams for three-particle 
processes (see Fig. 3), the positive definiteness of the amplitude squared is lo 

Results of the calculation of the structure function R with zeroth order plus 
0{as) expressions for F2 and Fi, as they appear in the expression for R are shown 
in Fig. 11. In the x > 0.25 region of interest the results typically lie a factor of 
6 above the experimental data for R. The absence of data for x > 0.8 in Fig. 11 
is due to the numerical instability encountered as Fi, appearing in a denominator 
in the expression for i?, goes to zero at large x. 

5. Summary 

In this work, we considered deep inelastic scattering of leptons off a proton 
in the statistical model proposed in ref. [2]. The interior of the nucleon is viewed 
as a thermalized assembly of up and down quarks and gluons. This enables one 
to incorporate features which are absent in the standard parton model. Quan- 
tum statistical correlations are incorporated through the use of Fermi-Dirac and 
Bose-Einstein distributions for initial-state quarks and gluons, respectively. Stim- 
ulated emission factors for final-state gluons and Pauli-blo eking factors for final- 
state quarks are incorporated. The propagation of particles through a many-body 
medium is taken into account by using thermal Feynman rules for propagators and 
vertices. The statistical model could also be seen as an attempt to describe the 
interior of the nucleon at a more fundamental level than that attained through the 
use of arbitrary parton distributions containing many parameters in the parton 
model. 

We chose to fit to a plot of F2 versus x, which is determined by parton 
distributions which reproduce a wide range of experimental data. Such parton 
distributions were published in 1990 by Kwiecinski et al. in [17]. They decide on 
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a parametrization of the parton distributions at Qq = —(^ = 4 GeV^ and evolve 
up in in order to test their parton distributions against experimental data. 

Our zeroth order in the strong coupling constant as calculation of the struc- 
ture function F2 versus x at — 4GeV^ produced a good fit to the data for 
X > 0.4. The values for temperature and chemical potential which give the best 
fit are T = 0.04 GeV and fiup = 0.21 GeV (with //down = Atup/2 3). Other com- 
binations of values for T and fi^p do not produce an improved fit to the data. In 
consequence, we investigate to what extent 0{as) corrections are able to extend 
the reproduction of the data to values of x smaller than 0.4. The correspond- 
ing zeroth order in as results for the structure function R are consistent with 
experimental data Fig. 7. 

The inclusion of the 0{as) corrections [7] in our numerical calculations allowed 
the extension of the reproduction of the data of Kwiecinski et al. to values of x 
smaller than 0.4. A fit could be obtained for values of a; > 0.25. This is compatible 
with our estimate for the limited range of validity of our theory {x > 0.21) due 
to a finite-size effect arising from the finite volume of the nucleon. The latter 
fit is possible for values of temperature and chemical potential in the immediate 
vicinity of T = 0.067 GeV and = 0.133 GeV only (with //down = /tup/23 and 
a, = 0.2). 

The latter values of parameters, however, do not fit the structure function 
i? = (Tl/(Jt to the experimental data presented in Ref. [18]. Even when taking 
into account that all measurements of R suffer from large experimental errors due 
to the weak dependence of the deep inelastic cross section for charged leptons on 
R, the size of the discrepancy remains unacceptable. This indicates a shortcoming 
of the statistical model in its present form to reproduce the structure function of 
the proton. 
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Appendix A : Final expressions for three- and four-particle processes 

In this appendix we summarize expressions derived for the phase space inte- 
grations of the matrix element squared from four- and three-particle processes. 

The way in which the final expressions for four-particle processes in eqs. (Bl) 
and (B2) and for three-particle processes in eqs. (CI) and (C2) contribute to the 
structure functions of the nucleon can be seen by considering eqs. (8), (9) and 
(15H18). 

In order to write the expressions in a simplified form, we introduce the ab- 
breviations 



A 


= {2K+po + v)/2 


a 


={2K + po-i^)/2 


B 


= {2K -po + iy)/2 


b 


={2K-pQ-v)/2 


C 


= {2K +po + qze{kok'oPo))/2 


c 


={2K -\-po - q:,e{kokQpo))/2 


D 


= {2K -po + q,e{kok'QPo))/2 


d 


={2K-pQ-q,e{kokM)/2 


E 


= {i^ + qze{kokQpo))/2 


e 


={u-q,e{kok'oPo))/2 (Al) 



where e is the sign function. 

In the calculation of the 0{as) corrections [9], all phase space integrations 
are analytically reduced to integrations over only the two energy variables ko and 
Po which are the energy components of the four-momentum assignments k and p, 
respectively, in Figs. 2 and 3. 

According to the mechanism by means of which collinear and infrared diver- 
gences cancel [9], it is appropriate to transform from the variables {ko,po) to the 
variables {K,po) where 

K^ko+'^ {A2) 

Appendix B : Final expressions for four-particle processes 

The convention concerning the sign of the energy of a particle discussed in 
Section 2, enables one to view all contributing four-particle processes simultane- 
ously in the two-dimensional energy plane shown in Fig. 4. The Latin letters 
in Fig. 4 serve to label the subregions over which one integrates in the K, pq 
plane. The numbers in Fig. 4 label the classes of subregions. The final analytical 
expressions are different for each class. The classes arise as follows: For some 
regions of the ko, /cq plane (/cq = ko + qo — po), the restrictions imposed by the 
energy-momentum conserving Dirac delta function in the phase space integration, 
causes one to replace the usual upper limit of 1 and/or the usual lower limit of 
— 1 for cos^ and/or cos^' (where 9 is the angle between q and k and 9' the angle 
between q and k') by expressions in terms of ko and /cq. The classes of subregions 
in the ko, kg plane are distinguished according to the combination of upper and 
lower limits for z and z' which are obtained for each class. In consequence, the 
integrations over z and z' will lead to different analytical results for each class. 
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The final results for four-particle processes after the cancellation of coUinear 
singularities are 
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with quantities as defined in Table 1 and 

1 e{ko) 



Ga 



{B2) 



Sa = 



2 2 

1 + e{po) 



(1 - 2nF{xk)) 
+ nB{\po\ 



^(l-2n44)) 



(53) 



These expressions are valid for all the four-particle processes with the K, po inte- 
gration region given in Fig. 4. 

The way in which these expressions contribute to the structure functions of 
the nucleon can be seen by considering eqs. (A2), (Bl) and (15)-(18). 

Appendix C : Final expressions for three-particle processes 

Similar to our discussions for the four-particle processes, both the quark 
and antiquark three-particle processes can be viewed simultaneously in the two- 
dimensional energy plane shown in Fig. 5. From the figure it can be seen that 
the K, pq integration region for three-particle processes is not subdivided into sub- 
regions as in the case of four-particle processes (see Fig. 4). For three-particle 
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processes, subdivisions of the energy plane need only to be considered when per- 
forming angular integrations in intermediate steps and these subdivisions disap- 
pear after the angular integrations. 

The final results for three-particle processes after the cancellation of collinear 
singularities are [Declaration of symbols appears after eq. (C2)]: 
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The corresponding final result for the longitudinal polarization of the virtual 
photon is: 
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+ S^F' { -e{B) 



2vA 



ipoql 
ln\BF' \ -e{B) 



2 

Z L 



ln\AF'\[q' + 4K'+pl-u']-e{B)^^^ 

\q^\ \ 1 a + A^^^^_^^^^ ^^1^^^ 



^°l4l 



Qz 



4 Ml 
12 



where 



1 <ko) 

2 2 



S3B = S2[l + 2nB{\po\)]^ 

SsF = [1 - 2nF{xk-p)] ^ 
SsF' = [1 - 2nF(a;fc/_p)] ^ 



;i - 2nF{xk)) 



^+ ^ ^ ^ (l-2nF(a;fc+,)) 



a^fc = |A;o| - //e(/co). 



Kf 



Dd 

{qz + v){2K -po - q^) 
{qz - v){2K -pQ + qz) 



Af' = 



Pi 



Dd 

{qz + v){2K -pQ + qz) 
{Qz - i^){2K -po - qz) 

Bf' — — 2" 
Po 



(C3) 

(C4) 

(C5) 

(C6) 
(C7) 

(C8) 

{C9) 

{CIO) 

{Cll) 
{C12) 

(CIS) 



and where, in the step functions $, 



y 



+ 2kou 
2koqz 



(CM) 



The factor S2 describes the statistical factors for the external either quarks or 
antiquarks in three-particle processes. The last term in each of eqs. (CI) and 
(C2) has been introduced [8] to make the integrations over positive and negative po 
separately convergent for the purpose of numerical calculation. From the definition 
of K it follows that K+po/2 is independent of po and therefore that the last term 
in each of eqs. (CI) and (C2) contributes zero when added at a positive and 
negative value of pq 
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Table 1 rdefinitions of the quantities appearing in eqs. (Bl) and (B2) [9] 



Class 


A4 


B4 


C4 




Fa 


Ga 


1 


1 


1 


1 


Po 


Po 


Po 


A{B-po) 


6(a + po) 


2 


B 


B 


Q 


B 


B 


B 






MD 


A(B-pq) 

J. J. y-*-^ j^yj J 


(d + e)(D + E) 


3 


7-) 
an 


LB 


La 


E 


E 




cD 




PoD 


A(c-E) 




4 


cB 


Be 


Poe 




c 


c 


ad 


VnE 


ad 


c 


A{c-E) 


(d -\- e)(a -\- Vn) 


5 


au 


Po-e 


au 


D 


n 


n 


BC 


c 


PoE 


{B-po){C-e) 


h{D + £;) 


6 


aB 


eB 


ea 


Po-e 


Po-e 




Cd 


PoC 


Pod 


{C-e){B-po) 


(a + po)(c^ + e) 


7 


a2 






a 


a 


a 


Cc 


4Cc 


2 
Po 


{C-e){c-E) 


6(a + ]3o) 
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Figure Captions 



Fig. 1 Processes contributing to deep inelastic scattering to first order in ag . 
Fig. 2 Diagrams for four-particle processes. 
Fig. 3 Diagrams for three-particle processes. 

Fig. 4 Regions of support for all four-particle processes in the K, po plane [9]. 

Fig. 5 Regions of support for all three-particle processes in the K,pq plane 

Fig. 6 Plots of F2 versus x at = 4 GeV^; solid line reproduces the parton dis- 
tributions of Kwiecinski et al. [17] fit) and dashed line represents the 
zeroth order in as expression given in eq. (10). 

Fig. 7 Plot of R versus a; at = 4 GeV^ for as = 0.2, T = 0.04 GeV, fXup = 0.21 

GeV and /idown = l^-up/'^^ 
Fig. 8 Plots of F2 versus a; at = 4 GeV^; short dashed line represents the zeroth 

order in expression given in eq. (10) and the long dashed line represents the 

zeroth order plus 0{as expression given in eqs. (8) and (15) with as = 0.2. 

The soHd fine is the 5_ fit of Ref. [17]. 

Fig. 9 Plots of F2 versus x a,t — 4 GeV^ as a function of temperature at fixed 
chemical potential fiup = 0.133 GeV according to the zeroth order plus 0{as) 
expression given in eqs. (8) and (15) with as = 0.2. The temperature is 0.043 
GeV for the lowest dashed curve 0.055 GeV for the next upper one, 0.067 
GeV for the third dashed line and 0.079 GeV for the highest dashed line. The 
soHd fine is the S_ fit of Ref. [17]. 

Fig. 10 Plots of F2 versus x at — 4 GeV^ as a function of temperature at fixed 
temperature T = 0.067 GeV according to the zeroth order plus 0{as) expres- 
sion given in eqs. (8) and (15) with as = 0.2. The chemical potential fj,up 
is 0.109 GeV for the lowest dashed curve 0.133 GeV for the next upper one, 
0.157 GeV for the third dashed line and 0.181 GeV for the highest dashed 
line. The solid line is the S_ fit of Ref. [17]. 

Fig. 11 Plot of i? versus x at ^ 4 GeV^ with as = 0.2, T = 0.067 GeV, Hup = 0.133 

GeV and Udown = i^up/'^^ 
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